Controlling paths of high-order harmonic generation by orthogonal two-color fields
Ma Ze-Hui1, 2, Zhang Cai-Ping1, 2, †, Ma Jun-Lin1, 2, Miao Xiang-Yang1, 2, ‡
College of Physics and Information Engineering, Shanxi Normal University, Linfen 041004, China
2Key Laboratory of Spectral Measurement and Analysis of Shanxi Province, Shanxi Normal University, Linfen 041004, China

 

† Corresponding author. E-mail: cpzhang880822@163.com sxxymiao@126.com

Project supported by the National Natural Science Foundation of China (Grant Nos. 11974229, 11404204, and 11947002), the Scientific and Technological Innovation Program of Higher Education Institutions in Shanxi Province, China (Grant No. 2019L0468), the Natural Science Foundation for Young Scientists of Shanxi Province, China (Grant No. 201901D211404), and the Innovation Project for Postgraduates of Shanxi Province, China (Grant No. 2019SY310).

Abstract

Controlling paths of high-order harmonic generation from is theoretically investigated by numerically solving the time-dependent Schrödinger equation based on the Born–Oppenheimer approximation in orthogonal two-color fields. Our simulations show that the change of harmonic emission paths is dependent on time-dependent distribution of electrons. Compared with one-dimensional linearly polarized long wavelength laser, multiple returns are suppressed and short paths are dominant in the process of harmonic emission by two-dimensional orthogonal ω/2ω laser fields. Furthermore, not only are multiple returns weaken, but also the harmonic emission varies from twice to once in an optical cycle by orthogonal ω/1.5ω laser fields. Combining the time–frequency distributions and the time-dependent electron wave packets probability density, the mechanism of controlling paths is further explained. As a result, a 68-as isolated attosecond pulse is obtained by superposing a proper range of the harmonics.

1. Introduction

Attosecond pulses[1,2] have attracted a lot of attention experimentally and theoretically, and it is widely used in studying electron dynamics and detecting the internal structure of atoms, molecules, and solids. Moreover, the attosecond pulses can be synthesized from high-order harmonics[3,4] which are emitted by the interaction of strong laser with atoms,[5,6] molecules,[79] and solids.[1013] The electron dynamics in harmonic generation from the atom and the small molecule with single electron are extensively investigated in one-dimensional (1D) laser field. But 1D linearly polarized light restricts the electron motion essentially to the light polarization plane where recollisions and harmonic emission occur nearly continuously over the laser pulse duration. The two-dimensional (2D) laser fields offer more possibilities for detecting the internal structure of molecules and studying the physical mechanism. At present, the 2D orthogonal fields have attracted a lot of attention.[1417] With orthogonally polarized two-color laser fields, Chen et al. researched temporal properties of high-order harmonic generation (HHG) from aligned molecules in orthogonal two-color laser fields[18] and probed the degrees of orientation of the nonlinear polyatomic molecules via odd and even harmonics.[19] He et al. studied HHG from in the orthogonal two-color pulses composed of strong fundamental frequency pulse and weak direct-current pulse, and found that the harmonics in the direction of weak direct-current pulse presented pure even harmonics, while those in the direction of strong fundamental frequency presented pure odd harmonics.[20] Gonzalez-ferez et al. proved that perpendicularly polarized two-color laser fields can be used to achieve stronger molecular orientation.[21] Lu et al. enabled the tomographic imaging of molecular orbital with single-shot measurement in the experiment by 2D processing of electron motion in the orthogonal fields.[22] 2D orthogonal fields can be used not only to detect molecular structure information and imaging information but also to control the paths of harmonic emission. The HHG is contributed by different trajectories, such as long and short paths, the interference between them is not conducive to synthesizing attosecond pulses. Many ways have been explored to synthesize attosecond pulses in the 1D linearly polarized laser field, such as few-cycle pulse[23] and the two-color scheme,[24] etc. Compared with the 1D laser field, the 2D orthogonal fields provide more possibilities to obtain the shorter and higher intensity attosecond pulses.

In this paper, we theoretically investigate the paths control of HHG from by numerically solving the Born–Oppenheimer (BO) time-dependent Schrödinger equation (TDSE). The results show that short paths are predominant while the long quantum paths are suppressed in the case of the 2D orthogonal ω/2ω fields. Moreover, the harmonic emission paths change from twice to once in an optical cycle (o.c.) in the case of the 2D orthogonal ω/1.5ω fields. This paper provides a simple and convenient method to control the harmonic emission paths, which can provide theoretical guidance for experiments.

2. Theoretical method

We assume that the fundamental field is linearly polarized along x axis, and the additional n-harmonic field along the y axis. With the molecular axis along the x axis, the 2D TDSE can be written as[2528]

Then, the Hamiltonian of the model molecule can be written as

(in atomic units of ). In the above expression, V(r) represents the Coulomb potential of the molecule which can be expressed as

with and in 2D cases. Here, R = 2 a.u. is the internuclear separation and Z is the effective nuclear charge. ζ = 0.5 is the smoothing parameter which is used to avoid the Coulomb singularity. We use a grid size of Lx × Ly = 400 a.u. × 250 a.u. with a grid spacing of Δx = Δy = 0.4 a.u. for the x and y axes, respectively. With the above parameters, the ground-state ionization potential of is 1.1 a.u. The electric field of 2D orthogonal fields can be written as

where Ex(t) = E0sin(ω0t) and Ey(t) = E1sin(0t). E0(E1) is the maximal amplitude of electric field Ex(t) (Ey(t)). ω0 is the laser frequency of Ex(t) lasting 5 o.c. Equation (1) is solved by the second-order split-operator method.[29] Once the wave function Ψ(t) is obtained, the dipole acceleration can be given by the Ehrenfest theorem[30]

In order to further understand the HHG, the time–frequency distribution is given by means of the wavelet transform[31]

3. Results and discussion

As we all know, the narrower attosecond pulse can be obtained by wider harmonic plateau. Compared with the short wavelength, the long wavelength is easier to extend the harmonic spectrum experimentally. It is the 1D linearly polarized long wavelength laser field (I = 5 × 1014 W/cm2, λ = 1200 nm) that is used in the paper, and corresponding the time–frequency distribution in Fig. 1. In the process of harmonic emission, there are not only first return, but also multiple returns as shown by white arrows. The underlying mechanisms of multiple returns have been explained in our previous works.[32,33] So, it is necessary to manipulate multiple returns, which make it difficult to synthesis the isolated attosecond pulse with the long wavelength pulse.

Fig. 1. Time–frequency distribution of the HHG corresponding to the case of I = 5 × 1014 W/cm2, λ = 1200 nm in 1D linearly polarized laser field.

Multiple returns will be eliminated effectively accompanied by the suppression of the long paths.[32] In essence, the motion time of electrons in the laser field dominates the long paths, which is mainly contributed by the early ionization and late return electrons. So the 2D orthogonal ω/2ω field (Ix = 5 × 1014 W/cm2, Iy = 5 × 1014 W/cm2, λx = 1200 nm, λy = 600 nm) in Fig. 2(a) is selected, which looks like a bowknot. Here A, A′, C, C′ (B, D) mean the ionization (return) moments of the electrons in one cycle, the red and blue arrows mean the trend of electric field.

Fig. 2. (a) The sketch of the electric field of the laser pulses with the combination of ω/2ω laser fields, Ix = 5 × 1014 W/cm2, Iy = 5 × 1014 W/cm2, λx = 1200 nm, λy = 600 nm. O is the origin of ExEy, A, A′, C, C′ (B, D) mean the ionization (return) moments of the electron, the arrows represent the trend of electric field; (b) the time–frequency distribution of the HHG corresponding to the case of panel (a); (c) and (d) the time-dependent electron wave packets probability density of x, y axes, respectively. The white lines mean that the excursion of the electron trajectories in the continuum is characterized by a distance x0 = E0/ω2.

Obviously, the multiple returns are suppressed in the time–frequency distribution as seen in Fig. 2(b). Besides, the harmonics are emitted twice in one cycle and mainly contributed by short paths. Figures 2(c) and 2(d) show the time-dependent electron wave packets probability density of x, y axes in 0 o.c.–2 o.c., respectively. To further explore the process of harmonic emission, the results in 0 o.c.–0.5 o.c. are taken as an example. There are twice ionization around point A′ (0.1 o.c.), A (0.4 o.c.), which are near the peak intensity of x, y axes in Fig. 2(a). When the electrons are ionized around point A′ (0.1 o.c.), they move farther and farther from the x axis and never come back, but return around point B′ (0.5 o.c.) in y axis as shown in Figs. 2(c) and 2(d). Furthermore, when the electrons are ionized around point A (0.4 o.c.), the electrons return to the parent nucleus simultaneously in the x, y directions around point B (0.8 o.c.) in Figs. 2(c) and 2(d). Only electrons returning to the parent nucleus simultaneously in the x, y directions, will contribute to HHG. This is similar to the other twice ionization around point C′ and C. Therefore, twice electron returns lead to twice harmonic emissions in one cycle.

To find the reason for suppression of multiple returns, classical computations are performed. The white lines in Fig. 2(c) mean that the excursion of the electron trajectories in the continuum is characterized by a distance x0 = E0/ω2 (82 a.u.). The short trajectories cannot go beyond x0 while the long trajectories always cross this value,[34] moreover, the HHG is also dominated by the instantaneous electric field direction parallel to the molecules axis.[22] As shown in the figures, the probabilities of electrons localized at the region between two white lines are overwhelming majority. It means that the long paths are weakened effectively in the process of HHG, which is consistent with Fig. 2(b).

To further optimize the harmonic emission, we choose the 2D orthogonal ω/1.5ω fields (Ix = 5 × 1014 W/cm2, Iy = 5 × 1014 W/cm2, λx = 1200 nm, λy = 800 nm) as displayed in Fig. 3(a), which looks like two-crossing heart. A, A′, B, B′, C, C′ (D, D′) indicate the ionization (return) moments of the electron, and arrows represent the trend of electric field. The time–frequency distributions of ω/1.5ω laser fields show that harmonic emission only once a cycle in Fig. 3(b). Figures 3(c) and 3(d) are the time-dependent electron wave packets probability densities in x, y axes in 0 o.c.–2.5 o.c. There are three ionization moments around point A (0.16 o.c.), B (0.5 o.c.), C (0.8 o.c.), respectively in one cycle. As discussed above, the harmonics are emitted once a cycle, corresponding to the process of C (0.8 o.c.)→ D(1.4 o.c.) in Figs. 3(c) and 3(d), respectively. We conclude that the harmonic emission in 2D orthogonal laser fields is strongly affected by time-dependent distribution of electrons. Furthermore, by superposing the harmonics from 120th order to 180th order, an isolated 68-as pulse is obtained as shown in Fig. 4.

Fig. 3. (a) The electric field of the laser pulses with the combination of ω/1.5ω laser fields in two cycles, Ix = 5 × 1014 W/cm2, Iy = 5 × 1014 W/cm2, λx = 1200 nm, λy = 800 nm, O is the origin of ExEy, A, A′, B, B′, C, C′ (D, D′) denote the ionization (return) moments of the electron, the arrows mean the trend of electric field; (b) the time–frequency distribution of the HHG corresponding to the case of panel (a); (c) and (d) the time-dependent electron wave packets probability density of x, y axes, respectively. The white line means that classical computations show that the excursion of the electron trajectories in the continuum is characterized by a distance x0 = E0/ω2.
Fig. 4. The time profiles of the attosecond pulses by superposing harmonics from the 120th order to the 180th order for Fig. 3(a).
4. Conclusion

In summary, we theoretically investigated the HHG process from by numerically solving the Born–Oppenheimer time-dependent Schrödinger equation. We present a scheme to control paths of the harmonic emission in the orthogonal two-color fields. The results show that the harmonic is emitted twice in one cycle by the 2D orthogonal ω/2ω laser fields, and it is mainly contributed by the short paths. Furthermore, harmonic emissions change from twice to once in one cycle by the 2D orthogonal ω/1.5ω laser fields. In contrast, HHG from 1D linearly polarized laser field is contributed by long paths, short paths, and multiple returns. The physical mechanism of the phenomenon is explained by the time–frequency analysis. In addition, by the time-dependent electron wave packets probability density, the law of motion of electrons in orthogonal fields is also investigated. The reason that the orthogonal two-color fields can choose the quantum paths and change the number of harmonic radiations is that the electron wave packet movement possesses different periodicities in the perpendicular direction. Moreover, an isolated 68-as pulse is directly obtained by superposing the proper range of the harmonics. Thus, our research can provide theoretical guidance for the experimental study.

Reference
[1] Corkum P B Krausz F 2007 Nat. Phys. 3 381
[2] Li P C Chu S I 2012 Phys. Rev. A 86 013411 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.86.013411
[3] Yuan K J Bandrauk A D 2019 Phys. Rev. A 100 033420 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.100.033420
[4] Milošević D B Becker W 2002 Phys. Rev. A 66 063417 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.66.063417
[5] L’Huillier A Schafer K J Kulander K C 1991 J. Phys. B: At. Mol. Opt. Phys. 24 3315
[6] Zhang J Qi T Pan X F Guo J Zhu K G Liu X S 2019 Chin. Phys. B 28 103204 https://cpb.iphy.ac.cn/EN/10.1088/1674-1056/ab4176
[7] Miao X Y Zhang C P 2014 Phys. Rev. A 89 033410 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.89.033410
[8] Li M Z Xu Y Jia G R Bian X B 2019 Phys. Rev. A 100 033410 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.100.033410
[9] Wustelt P Oppermann F Yue L Möller M Stöhlker T Lein M Gräfe S Paulus G G Sayler A M 2018 Phys. Rev. Lett. 121 073203 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.121.073203
[10] Yu C Jiang S C Lu R F 2019 Adv. Phys. X 4 1562982 https://www.tandfonline.com/doi/full/10.1080/23746149.2018.1562982
[11] Zuo R X Song X H Liu X W Yang S D Yang W F 2019 Chin. Phys. B 28 094208 https://cpb.iphy.ac.cn/EN/10.1088/1674-1056/ab3446
[12] Tao Z Chen C Szilvasi T Keller M Mavrikakis M Kapteyn H Murnane M 2016 Science 353 62
[13] Dejean N T Rubio A 2018 Sci. Adv. 4 eaao5207 https://advances.sciencemag.org/content/4/2/eaao5207
[14] Kitzler M Lezius M 2005 Phys. Rev. Lett. 95 253001 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.95.253001
[15] Zhang B Lein M 2019 Phys. Rev. A 100 043401 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.100.043401
[16] Milošević D B Becker W 2019 Phys. Rev. A 100 031401 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.100.031401
[17] Murakami M Korobkin O Horbatsc M 2013 Phys. Rev. A 88 063419 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.88.063419
[18] Chen C Ren D X Han X Yang S P Chen Y J 2018 Phys. Rev. A 98 063425 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.98.063425
[19] Yu S J Li W Y Li Y P Chen Y J 2017 Phys. Rev. A 96 013432 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.96.013432
[20] Xu L He F 2019 J. Opt. Soc. Am. B 36 840
[21] Mun J H Sakai H González-Férez R 2019 Phys. Rev. A 99 053424 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.99.053424
[22] Zhai C Y Zhang X F Zhu X S He L X Zhang Y F Wang B N Zhang Q B Lan P F Lu P X 2018 Opt. Express 3 2775
[23] Zhao S F Zhou X X Li P C Chen Z J 2008 Phys. Rev. A 78 063404 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.78.063404
[24] Du H C Wang H Q Hu B T 2010 Phys. Rev. A 81 063813 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.81.063813
[25] Lu R F Zhang P Y Han K L 2008 Phys. Rev. E 77 066701 https://journals.aps.org/pre/abstract/10.1103/PhysRevE.77.066701
[26] Miao X Y Zhang C P 2014 Laser Phys. Lett. 11 115301 https://iopscience.iop.org/article/10.1088/1612-2011/11/11/115301
[27] Hu J Han K L He G Z 2005 Phys. Rev. Lett. 95 123001 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.95.123001
[28] Miao X Y Du H N 2013 Phys. Rev. A 87 053403 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.87.053403
[29] Feit M D Fleck J A Jr Steiger A 1982 J. Comput. Phys 47 412
[30] Burnett K Reed V C Cooper J Knight P L 1992 Phys. Rev. A 45 3347
[31] Tong X M Chu S I 2000 Phys. Rev. A 61 021802 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.61.021802
[32] Zhang C P Pei Y N Jia X F Miao X Y 2017 Laser Phys. Lett. 14 015301 https://iopscience.iop.org/article/10.1088/1612-202X/14/1/015301
[33] Zhang C P Xia C L Jia X F Miao X Y 2016 Opt. Express 24 20297 https://www.osapublishing.org/oe/abstract.cfm?uri=oe-24-18-20297&origin=search
[34] Labeye M Risoud F Lévêque C Caillat J Maquet A Shaaran T Salières P Taïeb R 2019 Phys. Rev. A 99 013412 https://journals.aps.org/pra/abstract/10.1103/PhysRevA.99.013412